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Abstract 

We consider a new interacting particle system on the one- 
dimensional lattice that interpolates between TASEP and Toom's 
model: A particle cannot jump to the right if the neighboring site 
is occupied, and when jumping to the left it simply pushes all the 
neighbors that block its way. 

We prove that for flat and step initial conditions, the large time 
fluctuations of the height function of the associated growth model 
along any space-like path are described by the Airyi and Airy2 pro- 
cesses. This includes fluctuations of the height profile for a fixed time 
and fluctuations of a tagged particle's trajectory as special cases. 

1 Introduction 

We consider a model of interacting particle systems, which is a generalization 
of the TASEP (totally asymmetric simple exclusion process) and the Toom 
model. Besides the extension of some universality results to a new model, 
the main feature of this paper is the extension of the range of analysis to 
any "space-like" paths in space-time, whose extreme cases are fixed time and 
fixed particle (tagged particle problem), see below for details. 

Consider the system of particles Xi > ■ ■ ■ > xjy in Z that undergoes 
the following continuous time Markovian evolution: Each particle has two 
exponential clocks - one is responsible for its jumps to the left while the other 
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one is responsible for its jumps to the right. All 2N clocks are independent, 
and the rates of all left clocks are equal to L while the rates of all right clocks 
are equal to R. When the iih left clock rings, the ith particle jumps to the 
nearest vacant site on its left. When the ith right clock rings, the ith particle 
jumps to the right by one provided that the site Xj + 1 is empty; otherwise it 
stays put. The main goal of the paper is to study the asymptotic properties 
of this system when the number of particles and the evolution time become 
large. 

If L = then the dynamics is known under the name of Totally Asym- 
metric Simple Exclusion Process (TASEP), and if i? = the dynamics is 
a special case of Toom's model studied in [9] (see references therein too). 
Both systems belong to the Kardar-Parisi- Zhang (KPZ) universality class of 
growth models in 1 + 1 dimensions. 

Particle's jump to the nearest vacant spot on its left can be also viewed 
as the particle pushing all its left neighbors by one if they prevent it from 
jumping to the left. This point of view is often beneficial because it remains 
meaningful for infinite systems, and also the order of particles is not being 
changed. Because of this pushing effect we call our system the Pushing 
Asymmetric Simple Exclusion Process or PushASEP. 

Observe that for a A^-particle PushASEP with particles 
xi{t) > ■ ■ ■ > XN{t), the evolution of {xi,...^xm) for any M < is 
the M-particle PushASEP not influenced by the presence of the remaining 
N — M particles. This "triangularity property" seems to be a key feature of 
our model that allows our analysis to go through. 

Our results split in two groups - algebraic and analytic. 

Algebraically, we derive a determinantal formula for the distribution of 
the A^-particle PushASEP with an arbitrary fixed initial condition, and we 
also represent this distribution as a gap probability for a (possibly, signed) 
determinantal point process (see [12,16,17,21,22] for information on deter- 
minantal processes). The result is obtained in greater generality with jump 
rates L and R being both time and particle-dependent (Proposition [3]T]). The 
first part (the determinantal formula, see Proposition 12. ip is a generalization 
of similar results due to [2, 19, 20] obtained by the Bethe Ansatz techniques. 
Also, a closely related result have been obtained very recently in [10] using 
a version of the Robinson-Schensted-Knuth correspondence. 

Analytically, we use the above-mentioned determinantal process to study 
the large time behavior of the infinite-particle PushASEP with two initial 
conditions: 

1. Flat initial condition with particles occupying all even integers. 

2. Step initial condition with particles occupying all negative integers. 
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It is not obvious that the infinite-particle PushASEP started from these 
initial configurations is correctly defined, and some work needs to be done to 
prove the existence of the Markovian dynamics. However, we take a simpler 
path here and consider our infinite-particle system as a limit of growing 
finite-particle systems. It turns out that for the above initial conditions, 
the distribution of any finite number of particles at any finitely many time 
moments stabilizes as the total number of particles in the system becomes 
large enough. It is this limiting distribution that we analyze. 

We are able to control the asymptotic behavior of the joint distribution 
of Xn^{ti), . . .,Xnk{tk) with a;„,(0) > ■ ■ ■ > a;„^(0) and h > ■ ■ ■ >tk. It is the 
second main novel feature of the present paper (the first one being the model 
itself) that we can handle joint distributions of different particles at different 
time moments. As special cases we find distributions of several particles at a 
given time moment and distribution of one particle at several time moments 
(a.k.a. the tagged particle). 

In the growth model formulation of PushASEP (that we do not give 
here; it can be easily reconstructed from the growth models for TASEP and 
Toom's model described in [9] and references therein), this corresponds to 
joint distributions of values of the height function at a finite number of space- 
time points that lie on a space-like path; for that reason we use the term 
'space-like path' below. The two extreme space-like paths were described 
above - they correspond to ti = ■ ■ ■ = and rii = ■ ■ ■ = n^. 

The algebraic techniques of handling space-like paths are used in the 
subsequent paper [6] to analyze two different models, namely the polynuclear 
growth (PNG) model on a flat substrate and TASEP in discrete time with 
parallel update. 

Our main result states that large time fluctuations of the particle positions 
along any space-like path have exponents 1/3 and 2/3, and that the limiting 
process is the Airyi process for the flat initial condition and the Airy2 process 
for the step initial condition (see the review [11] and Section [23] below for 
the deflnition of these processes). 

In the PushASEP model, we have the fluctuation exponent 1/3 even in the 
case of zero drift. This is due to the asymmetry in the dynamical rules and it 
is consistent with the KPZ hypothesis. In fact, from KPZ we expect to have 
the 1/3 exponent when j"{p) ^ 0, where j(p) is the current of particles as a 
function of their density p, and j"(p) = —2{R + L/(l — p)^) for PushASEP. 

We flnd it remarkable that up to scaling factors, the fluctuations are 
independent of the space-like path we choose (this phenomenon was also 
observed in [7] for the polynuclear growth model (PNG) with step initial 
condition). It is natural to conjecture that this type of universality holds at 
least as broadly as KPZ-universality does. 
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Interestingly enough, so far it is unknown how to study the joint distribu- 
tion of Xniiti) and x„2(^2) with x„^(0) > a;„2(0) and ti < t2 (two points on a 
time-hke path); this question remains a major open problem of the subject. 

Previous results. For the TASEP and PNG models, large time fluctuation 
results have already been obtained in the following cases: For the step initial 
condition the Airy2 process has been shown to occur in the scaling limit for 
fixed time [14,15,18], and more recently for tagged particle [13]. For TASEP, 
the Airyi process occurs for fiat initial conditions in continuous time [4] and 
in discrete time with sequential update [Slwith generalization to the initial 
condition of one particle every d > 2 site^j. Also, a transition between the 
Airy2 and Airyi processes was obtained in [5]. These are fixed time results; 
the only previous result concerning general space-like paths is to be found 
in [7] in the context of the PNG model, where the Airy2 process was obtained 
as a limit for a directed percolation model. 

Outline. The paper is organized as follows. In Section [2] we describe 
the model and the results. In Proposition 12.11 the transition probability of 
the model is given. Then, we define what we mean by space-like paths, and 
formulate the scaling limit results; the definitions of the Airyi and Airy2 pro- 
cesses are recalled in Section fTM In Section [3] we state the general kernel for 
PushASEP (Proposition [m]) and then particularize it to step and flat initial 
conditions (Proposition 13.41 and 13. 6| ). In Section H] we first prove Proposi- 
tion 12.11 and then obtain the general kernel for a determinantal measure of 
a certain form (Theorem 14. 2p . which includes the one of PushASEP. Finally, 
the asymptotic analysis is the content of Section 
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2 The PushASEP model and limit results 



2.1 The PushASEP 

The model we consider is an extension of the well known totally asymmetric 
simple exclusion process (TASEP) on Z. The allowed configuration are like 
in the TASEP, i.e., configurations consist of particles on Z, with the con- 
straint that at each site can be occupied by at most one particle (exclusion 
constraint). We consider a dynamics in continuous time, where particles are 
allowed to jump to the right and to the left as follows. A particle jumps to 
its right-neighbor site with some rate, provided the site is empty (TASEP 
dynamics). To the left, a particle jump to its left-neighbor site with some 
rate and, if the site is already occupied by another particle, this is pushed to 
its left-neighbor and so on (push dynamics). 

To define precisely the jump rates, we need to introduce a few notations. 
Since the dynamics preserves the relative order of particles, we can associate 
to each particle a label. Let Xk{t) be the position of particle k at time t. We 
choose the right-left labeling, i.e., Xk(t) > Xk+i(t) for all A; G / C Z, t > 0. 
With this labeling, we consider > 0, k & I, and some smooth positive 
increasing functions a{t),b{t) with a(0) = 6(0) = 0. Then, the right jump 
rate of particle k is a{t)vk, while its left jump rate is h{t)/vk- 

In Proposition 12.11 we derive the expression of the transition probability 
from time t = to time t for N particles, proven in Section HI 

Proposition 2.1. Consider N particles with initial conditions a;i(0) = yi. 
Denote its transition probability until time t by 

G{xN, . . . , xi; t) = F{xi{t) = Xi,l<i< N\xi{0) = yul<i<N). (2.1) 

Then 

G{xN,...,xi;t) (2.2) 

N 



( n ^--^-e-^W'^-e-^W/""] det [Fk,i{xN+i-i - yN+i^k, ait), bit))] 



l<kA<N ' 



where 

F,,ix,a,b) = d^^-^ n£i(l-^A^+i~.^) ^..^a/.^ (2.3) 

^^^-^ro {[.^^il-VN+i-jz) 

where Fq is any anticlockwise oriented simple loop with including only the 
pole at z = 0. 
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2.2 Space-like paths 

The computation of the joint distribution of particle positions at a given 
time t can be obtained from Proposition 12.11 by adapting the method used 
in [4] for the TASEP. However, one of the main motivation for this work is to 
enlarge the spectrum of the situations which can be analyzed to what we call 
space-like paths. In this context, space-like paths are sequences of particle 
numbers and times in the ensemble 

S = {{nk,tk), k > l\{nk,tk) -< (nfc+i,4+i)}, (2.4) 

where, by definition, 

ijii^ti) -< {nj,tj) if Tij > ni,tj < ti, and the two couples are not identical. 

(2.5) 

The two extreme cases are (1) fixed time, tk = t for all k, and (2) fixed 
particle number, Uk = n for all k. This last situation is known as tagged 
particle problem. Since the analysis is of the same degree of difficulty for any 
space-like path, we will consider the general situation. 

Consider any smooth function vr, = 7r{w^), in the forward light cone 
of the origin that satisfies 

k'l < 1, \w^\<'k{w^). (2.6) 

These are space-like paths in R x ]R_|_, see Figure [TJ The first condition (the 
space-like property) is related to the applicability of our result to sequences 
of particles in S. The second condition just refiects the choice of having t > 
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and n > 0. Time and particle number are connected with the variables 
and by a rotation of 45 degrees. To avoid unnecessary \/2's, we set 



(2.7) 



We want to study the joint distributions of particle positions in the limit 
of large time, where universal processes arise. Since we consider several 
times, we can not simply use t as large parameter. Instead, we consider a 
large parameter T. Particle numbers and times under investigation will have 
a leading term proportional to T. In the {w^,w^) plane, we consider 
around OT for a fixed 6, while = Tn{w^/T). From KPZ we know that 
correlations are on T^/^ scale. Therefore, we set the scaling as 

' w\u) = er - uT^/^ 

«;0(m) = n{e)T - n'{e)uT^/^ + \Ti"{e)u''T^I\ ^ ' ' 

Notice that w^{u) is equal to Ttt{w^{u)/T) up to terms that remain bounded, 
and they become irrelevant in the large T limit, since the fluctuations grow 
as T^/^. Coming back to the {n,t) variables, we have 

t{u) = {7v{e) + e)T - {7i'{e) + i)uT^/'' + l-n" {e)u^T^'\ 

n{u) = [{'R{e)-e)T + {i-Tx\e))uT^'^ + \'K"{e)u^T^/'']. (2.9) 

In particular, setting ti{9) = 1 — 6 we get the fixed time case with t = T, 
while setting 7i{9) = a + 9 we get the tagged particle situation with particle 
number n = aT. 



2.3 Scaling limits 

Universality occurs in the large T limit. In Proposition 13.11 we obtain an 
expression for the joint distribution in the general setting. For the asymptotic 
analysis we consider the case where all particles have the same jump rates, 
i.e., we set 

Vk = l for all ke I. (2.10) 

Moreover, we consider time-homogeneous case, i.e., we set a{t) = Rt and 
b{t) = Lt for some R,L>0 (for time non-homogeneous case, one would just 
replace R and L by some time-dependent functions). Two important initial 
conditions are 

(a) flat initial condition: particles start from 2Z, 

(b) step initial condition: particles start from Z_ = {. . . , —3, —2, —1}. 
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In the first case, the macroscopic limit shape is flat, while in the second case 
it is curved, see [11] for a review on universality in the TASEP. For TASEP 
with step initial conditions and particle-dependent rates Vk, the study of 
tagged particle has been carried out in [13]. 



Flat initial conditions 

For the fiat initial condition, it is not very difficult to get the proper scaling 
limit as T — > oo. The initial position of particle n{u) is —2n{u) and during 
time t{u) it will have traveled around vt(u), where v is the mean speed of 
particles, given by 

V = -2L + R/2. (2.11) 

The reason is that the density of particle is 1/2 and the particles jumps 
to the right with rate R but the site on its right has a 1/2 chance to be 
empty. Moreover, particles move (and push) to the left with rate L but 
typically every second move to the left is due to a push from another particle. 
Therefore, the rescaled process is given by 

u ^ Xt[u) = -y^^ , (2.12) 

where n{u) and t{u) are defined in (12.91) . The rescaled process Xt has a 
limit for large T given in terms of the Airyi process, Ai (see [4,5, 11] and 
Section for details on ^i). 

Theorem 2.2 (Convergence to the Airyi process). Let us set the vertical 
and horizontal rescaling 

o _ ..o r , pw^..^ , .^^l/3 o _ 4((8L + i^)(7r(g) + g))^/3 

(2.13) 

Then 

lim Xriu) = S^Ai{u/Sh) (2.14) 

T — *oo 

in the sense of finite dimensional distributions. 

The proof of this theorem is in Section The specialization for fixed 
time t = T is 

S^ = {8L + Rr\ 5. = ^^^^, (2.15) 

and the one for tagged particle n = aT at times t(u) = T — 2uT^/^, obtained 
by setting 6* = (1 — a)/2, is 

S. = i8L + R)y\ 5. = ^^^^. (2.16) 
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Figure 2: Parametric plot of for L = 1, i? = 4. 



Step initial condition 

The proper rescaled process for step initial condition is quite intricate. De- 
note by jSt the typical position of particle with number around at at time 
t. In the situations previously studied in the literature, there was a nice 
function /? = /?(«). In the present situation this is not anymore true, but we 
can still describe the limit shape. More precisely, a and f3 are parametrized 
by a /i G (0, 1) via 



a{fi) = il-fi)\R + L/n^), (3{fi) = -{{l-2fi)R + L/fi^). (2.17) 



The parameter /i comes from the asymptotic analysis in Section [5121 where it 
represents the position of the double critical point. To see that it is a proper 
parametrization, we have to verify that for a given point on the space-like 
curve {9, 7i{9)) there corresponds exactly one value of /i. From (12. 9p we have 
n ~ t{7i{6) — 9)/{'k{9) + 9) and, since we have set n ~ at, we have 



For any given 9, there exists only one /i such that (12.181) holds, because a is 
strictly monotone in /x. Some computations are needed, but finally we get 
the rescaling of the position function of m, namely. 



a{n) 



7i{9)-9 
7i{9) + 9' 





(2.19) 
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where 



Ox ■■ 



1 + in'iO) + 1) (^i? - + (1 - (2-20) 

, (vr'(^)(l-a(^))-(l + «(^)))2 



c^2 = i7r"(0) + 



L ]_ 



4(l-/x)3(7r(^) + 0)(i? + L//i3)' 



The rescaled process is then given by 



u ^ = , (2.21, 

with n{u) and t{u) given in 02.91) . Define the constants 

(7r(^) + e)(/? + L//i3) 



/i(l - /i) 



^ ^ (7r-(g) + l)(i? + L/^^) 7r-(g)-l 
^ 2/i 2/i(l-/i)2' ^ ■ ' 

Then, a detailed asymptotic analysis would lead to, 

lim Xt{u) = ^1^1^^ A2{niK,o'^^^u) , (2.23) 

in the sense of finite dimensional distributions, where A2 is the Airy2 process 
(see [11,14,18] and Section [231 for details on ^2)- As for the fiat PNG, special 
cases are tagged particle and fixed time. In Section 15.21 we obtain fl2.23p 
by looking at the contribution coming from the series expansion around a 
double critical point. To get (12.231) rigorously, one has to control (1) the 
error terms in the convergence on bounded sets and (2) get some bounds to 
get convergence of the Fredholm determinants. This is what we actually do 
in the fiat initial condition setting. 



2.4 Limit processes 

For completeness, we shortly recall the definitions of the limit processes Ai 
and A2 appearing above. The notation Ai(a;) below stands for the classical 
Airy function [1]. 

Definition 2.3 (The Airyi process). The Airyi process ^1 is the process 
with m-point joint distributions at Ui < U2 < ■ ■ ■ < Um given by the Fredholm 
determinant 

m 

P( []{Ai{uk) < Sk}) = det(l - XsKA,Xs)LH{uu...,ur^}>cu), (2.24) 

k=l 
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where Xs{uk,x) = l{x > Sk) and the kernel K^^ is given by 

Kai{ui,si;u2,S2) = , ^ =exp -— ^ ^— l{u2 > ui) 

V47r(M2-Mi) V 4(u2-Mi)y 

+Ai(si + S2 + {u2 - uif) exp ( {u2 - mi)(si + S2) + ^(^2 - ui)A . (2.25) 



Definition 2.4 (The Airy2 process). The Airy2 process A2 is the process 
with m-point joint distributions at ui < U2 < . . . < Um given by the Fredholm 
determinant 

m 

P( f]{A2{uk) < Sk}) = det(l - XsKA,Xs)L\{uu...,ur^}xu), (2.26) 

k=l 

where Xs{ukiX) = l{x > Sk) and the kernel Kjs,^ is given by 

/^^e-^(---)Ai(.i + A)Ai(.2 + A), U2>u,, 

Ka.Au.,S,,U2,S2)-{ .J^ ,-A(„.)Ai(., + A)Ai(.2 + A), U2<U,. 

(2.27) 



3 Finite time kernel 

In this section we first derive an expression for the joint distributions of par- 
ticle positions in a finite system. They are given by Fredholm determinants 
of a kernel, which is first stated for general jump rates and initial positions. 
After that, we specialize to the cases of uniform jump rates in the case of step 
and fiat initial conditions. Flat initial conditions are obtained via a limit of 
finite systems. 

3.1 General kernel for PushASEP 

To state the following result, proven in Section HI we introduce a space 
of functions Vn- Consider the set of numbers and let 

{ui < U2 < . . . < Uy] be their different values, with being the multiplicity 
of Uk {vk is the jump rate of particle with label k). Then we define the space 

K = span{x'u^, 1 < A; < z/, < / < ftfc - 1}. (3.1) 

Recall that the evolution of particle indexed by n is independent of the 
particles with index m > n. 
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Proposition 3.1. Consider a system of particles with indices n = 1,2,... 
starting from positions yi > 1/2 > ■ ■ ■■ Denote by the position of particle 
with index n at time t. Then the joint distribution of particle positions is 
given by the Fredholm determinant 

m 

P{r\Mt,) > s,}) = det (1 - X.i^X.)..({(„„M,...,(n.,t.)}x^) (3-2) 
fc=i 

with ((ni,ti), . . . , {nm,tm)) e S, and Xs{{nk,tk)){x) = l{x < Sk). The kernel 
K is given by 

k=l 

(3.3) 

where 

vi/;^^^(a;) = — / d^^-^'-^e"W/-+^W- ~ ""'l' , / = 1,2, 



27ri /p^ " " " (1 - viz) ■ ■ • (1 - viz) ' 

(3.4) 

the functions are uniquely determined by the orthogonality rela- 

tions 

E '^n-i{^)^n-k{^) = ^KU l<kj<n, (3.5) 
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and by the requirement span{$"lj(x), I = 1, . . . ,n} = Vn- The first term 
(IX^j is given by 

- 2ni 1 zy-^-^^ (1 - .„,+,.) ... (1 - .„,.) 

(3.6) 

r/ie notation Tq stands for any anticlockwise oriented simple loop including 
only the pole at 0. 

3.2 Kernel for step initial condition 

We set all the jump rates to 1: Vi = V2 = ■ ■ ■ = 1- The transition function 
03.61) does not depend on initial conditions. It is useful to rewrite it in a 
slightly different form. 

Lemma 3.2. The transition function can be rewritten as 

0((ni,ti),(n2,t2))(^^^) (3.7) 



oa{t2)w+b{t2 ) /w ^ [("1 '*i )^ ('^2 ,t2)] • 
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Proof of Lemma \3.2[ The proof follows by the change of variable z = 1/w in 



Lemma 3.3. Let yi = > 1. Then, the functions $ and \& are given by 

<|.;'*(x) = — ^ dz e-W-^W/-. (3.8) 

^ 27ri /r^ (2; - 1)^+^ 

Proof of Lemma \3.3i \E'^'*(x) comes from the change of variable z = 1/w in 
(13.41) ■ For /c > 0, the pole at w = 1 is irrelevant, but in the kernel enters 
also for negative values of k. 

We have to verify that the function satisfy the orthogonal con- 

dition (13. 5p and span the space Vn given in (13.11) . For vi = ■ ■ ■ = Vn = I, 
Vn = span(l, X, . . . , x"~^). By the residue's theorem, the function <l>"'*(x) is 
a polynomial of degree j in x. Thus, span($J'*(x), j = 0, . . . , n — 1) = Vn- 

The second step is to compute Xlxez '^j'*(^)^fc'*(^) for < j. A; < n — 1. 
We divide it into the sum over x > and the one over x < 0. We have 



x>0 x>0 ^ ' "^^1 "'^o 



a{t)w+b{t)/w j^\fc ^x+n 



^a{t)z+b(t)/z _ iy+1 yjX+n+1 ' 

(3.9) 



We choose the paths Fq and Fi satisfying \z\ < \w\, so that we can take the 
sum inside the integrals and use 

^X 1 

E^ = ^' (3-10) 



W^^^ W — Z 
x>0 

to get 

If f ^a{t)w+h{t)/w _ ]^ 

= ——^* dzi <^w———r^r — T-^ , (3.11) 

where the subscript z in Fo,^ reminds that z is a pole for the integral over w. 
Next consider the sum over x < 0. We have 



$:$;''(x)vi/f(x)=$:-i-/ d^/ d. 

x<0 x<0 ^ ' "^^0 "Ti 



a{t)w+b{t)/w -^\k ^x+n 



Qa{t)z+b{t)/z _ iy+1 ^x+n+1 ' 

(3.12) 
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This time we choose the paths Tq and Fi satisfying \z\ > \w\ and then take 
the sum inside the integrals. Using 

E^ = — (3-13) 

x<0 

we obtain 

where now w is a pole for the integral over z. Thus, 

(a;) = dsnD + ra. (3.15) 

We can deform the paths of integration in fl3.14l) so that they become as the 
integration paths of (13.111) up to correcting the contribution of the residue 
dX z = w. Thus we finally get, for < j, < n — 1, 

^f{x)n\^) = ^J - If-^-' = 5,,,. (3.16) 

As a side remark, the same computations would not hold for k,j < 0, for 
which $"'*(x) = 0, because it is not possible to choose the paths with \z\ > \w\ 
without introducing an extra pole at 2; = 0. □ 

Proposition 3.4 (Step initial condition, finite time kernel). 

The kernel for Hi = —i, i>l, is given by 

K{{ni, ti),xi; (na, ta), 2:2) (3.17) 

If 1 f W ^a{ti)w+b{ti)/w 



^ (27ri)2 /p^ '^^/n e^(*2)/^+'^(*2)^ ^y^i+m+i (1 _ ^)n2 w - z' 

The contours Fq and Fi include the poles w = Q and z = 1 and no other 
poles. This means in particular that Fq and Fi are disjoints, because of the 
term l/{w — z). 



Proof of Proposition \3.4\ Consider the main term of the kernel, namely 

I 1 1 1 ^ To i 



k=l k=l ''^0.1 

X 



1 r X2+n2 

_L (/ dz-—- ^e-"(*2)"-*(*2)/^ (3.18) 
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First we extend the sum to +00, since the second term is identically equal 
to zero for k > n2. We choose the integration paths so that l^; — 1| < |w — 1|. 
Then, we can take the sum inside the integral. The fc-dependent terms are 

J^iE^.^. (3,19, 

k>l ^ ' 

Thus, we get 

1 r r ^a(ti)w+b{ti)/w _ j^'jni ^X2+n2 ^ 

(27ri)2 /p^ ^ ea{t2)z+b{t2)/z ^xi+m+l _ X)n2 w - z' 

(3.20) 

Notice now we have a new pole aX w = z, but the one at w = 1 vanished. 
The contribution of the pole at w = 2; is exactly equal to the contribution 
of the pole at z = 1 in the transition function (13. 7p . Therefore in the final 
result the first term coming from (13.71) has the integral only around 2; = 0, 
and the second term is (I3.20p but with the integral over w only around the 
pole at tt; = and does not contain z. Finally, a conjugation by a factor 
(— l)"i""2 gives the final result. □ 



(EIHl) = 



3.3 Kernel for flat initial condition 

The kernel for the flat initial condition is obtained as a limit of those for 
systems with finitely many particles as follows. We first compute the kernel 
for a finite number of particles starting from yi = —2i, i > 1. Then we shift 
the focus by particles, i.e., we consider particles with numbers N + rii 
instead of those with numbers nj. For any finite time t, we then take the 
N —>■ 00 limit, in which the deviations due to the finite number of particles 
on the right tend to zero. The limiting kernel is what we call the kernel for 
the flat initial condition {yi = —2i with i e Z). 

For this case we also consider the homogeneous jump rates, Vi = V2 = 
■ ■■ = 1. 

Lemma 3.5. Let y^ = —2i, i > 1. Then, the functions $ and "if are given 
by 



^"'*(a;) = (f> dw — — 



(W(W - l))\ ^(,)^^fe(,)/, 

ro.i 



'J 



{2z — '^)z''^'^'^ a{t)z~b(i)/z 



2m /r, {z{z - l))^+i 



15 



Proof of Lemma \3.5[ The proof is almost identical to the one of Lemma I3.3[ 
The only difference is that contribution of the residue is in this case given by 

J2^-'\x)^r{x) = ^J dz{2z - l){z{z - = 5,,, (3.22) 

by the change of variable w = z{z — 1). □ 

Proposition 3.6 (Flat initial conditions, finite time kernel). 

The kernel for Hi = —2i, i^TL, is given by 

K{{ni,ti),xi; in2,t2),X2) 

I r I f W ^a{ti)w+b(ti)/w 



— 1 f ga(ti)(l-2:)+6(ti)/(l-z) ^n\+n2+X2 
+ TTIT ,.r,(i^\-.^h(i.My T^ .X„,4-r,o^^.,-H • (3.23) 



Proof of Proposition \3.(A The first step is to get the kernel for yi = —2i, 
i > 1. This step is similar to the one of Proposition 13.41 The integration 
paths are taken to satisfy \z{z — 1)\ < \w{w — Then this time, the sum 
in k is 

^ (z(z - 1))>'-^ 1 

{w{w - 1))'' {w-z){w-l + z) ^' ' 

and we get 

]_ /" r ^a(ti)w+b{ti)/w _ j^Nni ^X2+n2 2z — 1 

dz (h dw- 



(27ri)2 ^ ea(t2)2+fe(t2)A (^ - l)n2 yjxi+m+i _ _ 1 _^ ^) • 

(3.25) 

Notice that the pole for w = 1 is now replaced by two simple poles, one at 
w = z and one at w = 1 — 2;. The pole at w = z cancels with the one at 
2; = 1 of (13. 7p . Therefore, the kernel for y^ = —2i, i > 1, is given by 

~ 2^i ) e^{t2)w+bit2)/u,hiniM)<n2,t2)] 

Qa{ti)w+b(ti)/w _ j^^ni ^X2+n2 2z — 1 



dz (b dw- 



(27ri)2 Jp^ Jp^ (,a{t2)z+b(t2)/z (^^ _ ]^)n2 ^x'l+m+l (^^y _ ^j) _ 1 _j_ ' 

(3.26) 

At this point, we pick a large N and shift the focus to particles around 
the A^th one. Accordingly, we shift the positions by —2N. More precisely, in 
fl3.26p . we replace 

ni^ni + N, Xi- 2N. (3.27) 
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Then we get the kernel K = Kq + Ki + K^^'> with ((rii, ti), xi; {n2,t2), X2)- 
entries given by 

If I f W ^a{ti)w+b(t^)/w 



Iff ^aiti)w+b{ti)/w _ ]^\ni+Af _^x2+n2-Ar 2z - 1 

i^''' = ^ Q) dz (h dw- 



(27ri)2 Jr,, e'^(*2)^+&(t2)A _ i)n2+Jv ^zi+m-Tv+i _ _ i + 

The terms Kq and are independent of N, while i^*^^) is not. We need 
to show that in the — > oo limit, the contribution of K^^^ vanishes, in 
the sense that the Fredholm determinant giving the joint distributions of 
Proposition 13. II converges to the one with kernel Kq + Ki. 

The Fredholm determinant fl3.2l) is projected onto Xj < Sj. Therefore, for 
any given si, . . . ,Sm we need to get bounds on the kernel for Xj's bounded 
from above, say for Xi < £ for an £ e Z fixed. In the simplest case of 
pure TASEP dynamics {b{t) = 0), the limit turns out to be easy because 
for Xi + 111 < N, the pole at = vanishes. However, in our model, 
b{t) is generically non-zero and the integrand has an essential singularity 
at w = 0. In what follows, we choose the indices k G {1, . . . yfn} so that 
{nk,tk) -< (n/c+i, tfc+i), k = 1, . . . ,m — 1. Also, we simplify the notation by 
writing k G {l,...,m} instead of {nk,tk) in the arguments of the kernel. 
Then, the Fredholm determinant becomes 

(_l)n 

l<k,l<n ■ 

n>0 ii,...,i„=l xi<Si-^ Xn<Si^ 

(3.29) 

We apply the following conjugation of the kernel, which keeps unchanged the 
above expression, 

K{ik, Xfc; lu xi) = K{ik, Xk, lu a;,)e-'^'— '=^'=e(^'-^'=)/2. (3.30) 

Using the bound of Lemma 13.71 for any choice of e in (0, (8m)~^] and for 
Xk^xi bounded from above, we have 

\Ko{ik,Xk]ihXi)\ < const e^""', 

\ki{ik,Xk]iuXi)\ < const e(^'+^'=)/*^ < const e^^S (3.31) 
\K^^\ik,Xk;iuXi)\ < const e("'+"'=)/®fi:^ < const e""'/€^, 

where k G [0, 1). In the above bounds, we use the same symbol 'const ' for 
all the constants. With the choice of ordering of the {nk,tk)^s, we have that 
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Kq = = Kq if ii < ik, thus the bound holds trivially. For the case ii > ik, 
Lemma 13.71 implies the estimate 

l^ol < const e"^'=(^'-^'=)e-l^'-^'=l/V^S (3.32) 

for Xk,xi bounded from above. The bound in (13.311) is then obtained by 
choosing e < (4m)~^, since then eik < sm < 1/4. The other bounds on Ki 
and K^'^'^ are satisfied for em < 1/8. 

Therefore, the summand in the multiple sums of (13.291) is uniformly 
bounded by 



-1)" 

det [K{ik,Xk;ii,xi] 



n\ \ / l<kl<n 



< i-e=(^i+-+^")const"(l + /€^)"n"/^ 
n\ 



(3.33) 

the term n"/^ being Hadamard bound on the value of a n x n determinant 
whose entries have modulus bounded by 1. Since k < 1, replacing 1 + 
by 2 yields a uniform bound, which is summable. Thus, by dominated con- 
vergence we can take the N —>■ oo limit inside the Fredholm series. Since 
K < 1, we have limTv^oo K^'^^ = 0, thus the result is proven. Finally, just 
for convenience, we conjugate the kernel by (— l)"^^'^^^ which however has no 
impact on the Fredholm determinant in question. □ 

Lemma 3.7. Let Kq, Ki, K^^^ be as in lis. 28\) . Then, for Xi,X2 < i, we 
have the following bounds. 

\Ko{{n^,t,),x^■, (n2,t2),X2)| < const e(^-^'^)/2g-|x.-.i|/4j^^^^^^^^^^^^_^^^^^ 

\Ki{{n,,t,),xi;{n2,t2),X2)\ < const e^"'^/^e^''+'^'^/\ (3.34) 
\K^''\{ni,ti),Xi;{n2,t2),X2)\ < const e("i-"2)/2g(xi+x2)/4^Ar^ 

for some k G [0, 1). The constants const and k are uniform in N and depend 
only on i and Ui, ti 's. 

Proof of Lemma \37i\ For Kq and X2 — xi > 0, we can just choose the integra- 
tion path as Tq = {\w\ = e~^}, from which we have \Kq\ < const e"*-^'^"^^-* < 
const e~(^2-a;i)3/4_ ^Yie case xi — X2 > 0, we choose the integration path as 
Fo = {\w\ = e-i/^}. Then, \Kq\ < const e(^i-^2)/4_ 
For Ki, we choose Fi = {|1 — 2;| = e~^}. Then, 



Ki\ < const — ' , . (3.35) 



maxpi 


z 


X2 


minpj 


1 




z 


XI 



Along Fi, |1 — 2;| is constant, thus (minr^ |1 — 2;|^'^) ^ = e^^^ < const e^^^/^ for 
Xi bounded from above. Remark that const depends on the upper bound. 
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i, for Xl- In this case, we can take const = e^^/^. Also, for X2 bounded from 
above, maxp^ \z\^^ < const (1 — l/e^Y^ < const e"^'^/^. 

For K^^\ we use the path Tq = {\w\ = e^^} and Fi = {|1 — z\ = e~^}. 
As required, these paths do not intersect because 1/e^ < 1 — 1/e^. Then, 

|ir(^)|<const"^""^^'"f^^, K="^""^° \w{w^l)\_ 

For Xl, X2 bounded from above, we have maxr^ \z\^^ < const e~^^/^ (as above), 
and (minpg |if|^'^)~^ = e^^^ < const e^^'^^^. Finally, it is not difficult to obtain 
K = (1 + l/e^)/(l — e^) = 0.159. . ., since the maximum of \w{w — 1)| is 
obtained at w = —e~^ and the minimum of \z{z — 1)| at 2; = 1 — e~^. □ 



4 Determinantal measures 

In this section we first prove Proposition 12. 1[ Then, we use it to extend 
the measure to space-like paths. More precisely, we first obtain a general 
determinantal formula in Theorem 14. 1[ Then, in Theorem 14.21 we prove that 
the measure has determinantal correlations and obtain an expression of the 
associated kernel. 

Proof of Proposition l2J[ We first prove that the initial condition is satisfied. 
We have 

F.Kx,0) = ^ / d..-^nSl^^^^^. (4.1) 

(a) Fk,i{x, 0) = for a; > 1 because the pole at 2; = vanishes. 

(b) Fk^i{x, 0) = for A; > Z and x < I — k, because then 

Fk,i{x, 0) = / dzz^-\l - Viz) Vk-iz) (4.2) 



27ri 



ro 



and the residue at infinity equals to zero for x < I — k. 

Assume that x^ < ■ ■ ■ < Xi. If x^v > y^, also > t/at for / = 1, . . . , A^— 1. 
Thus Fi^i{xiy+i^i — niy, 0) = using (a). Therefore G{xiy, . . . , Xi; 0) = 0. On 
the other hand, if xat < i/n, then xjy < i/k — N + k, k = 1, . . . , N — 1. Thus 
Fk,i{xN — VN+i-k, 0) = using (b) and the fact that xn — i/N+i-k < I — k. 
Therefore we conclude that G{xn, . . . , xi; 0) = if xat 7^ i/n- For xn = i/n, 
Fi,i(0,0) = 1 and by (a) F^ix^+i-i - y^v, 0) = for / = 2, . . . , AT. This 
means that 

G{xn, . . . , xi; 0) = ^^^.^^^(xAr.i, . . . , xi; 0). (4.3) 
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By iterating the procedure we obtain 

G{xn, . . . ,xuO) = Y\_S, 



N 



^k,yk- 



(4.4) 



k=l 



Notice that the prefactor in (12.21) is equal to one at t = 0. 

The initial condition being settled, we need to prove that (12. 2p satisfies 
the PushASEP dynamics. For that purpose, let us first compute 



dFk,iix,t) 
dt 



a{t)Fk,i{x - 1, t) + h{t)Fk,i{x + 1, t). 



(4.5) 



from which it follows, by differentiating the prefactor and the determinant 
column by column. 



dG{xN, . ..,xi]t) 
dt 



-(^a{t)^Vk + hit) ^ —^G{xn, ...,xi;t) 
k=i k=i ^'^ 

N 

+a{t)^VkG{...,Xk-l,...;t) (4.6) 

k=l 
N 

+6(t)^-G(...,a;, + l,...;t). 



1=1 



To proceed, we need an identity. Using 



+ z 



1 - Vn+1-IZ 1 - Vn+1-IZ 

it follows that 

Ffc,z+i(x,t) = Fk,i{x,t) + VN+i-iFk,i+i{x + l,t). 
Therefore, for j = 2, . . . , iV, by setting yk = i/N+i-k, 



(4.7) 



(4.^ 



— det 



N+l-l Fk,l{^N+l-l — Vk, t) 

"^j' Fk,N+i-j{xj - yk, t) v'^'_^Fk,N+2-j{xj-i - ijk, t) ■ 



l<k,l<N 



(4.9) 



Here does not depend on the x/s. 
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Using fl4.8l) we have 

Vj'_iFk,N+2-j{xj - Vk, t) (4.10) 



^^jii^yfc,iv+i-i(a;j -yk,t)+ v''j^X^Fk^M+2-j{.Xj + l-yk,t) 



Using this identity in the previous formula, the first term cancels being pro- 
portional to its left column, and the second term yields 

G{. . . , Xj,Xj-i = Xj, . . . ; t) = — —G{. . . , Xj, Xj^i = + 1, . . . ; t). (4.11) 

Vj-i 

With fl4.1ip we can go back to fl4.6p . First, consider all the terms in (14. 6 p 
which are proportional to a{t). They are given by 

N N 



- J2 ^kG{. ..■,t) + J2 ^kG{. . . , Xfc - 1, . . . ; t) (4.12) 

k=l k=l 

N 

= ^y^Gi...;t)-Y,Vk{l-6,,_,,,,^,)Gi...-t) (4.13) 

k=2 

N-1 

+vnG{xn - 1, . . . ; t) + ^ ffe(l - 4,+i,xJG(. . . , Xfc - 1, . . . ; t) (4.14) 

k=l 

N 

- ^ VkG{. ..,Xk, Xk-i = Xk + I, . . .;t) (4.15) 

k=2 
N~l 

+ ^ VkG{. . . , Xk+1 = Xk, Xk,...]t). (4.16) 

fc=i 

The notation means that the first term of f l4.12p has been subdivided into 
fl4.13p . which contains non-zero terms when Xk-i 7^ Xk + 1, and 04.151) . whose 
terms are non-zero only when Xk-i = + 1. Similarly for the second term 
of (14.121) . By using (14.111) and shifting the summation index by one, we get 
that (14.161) equals 

N 

y^Vk-iG{. . .,Xk,Xk-i =Xk + l,.. .;t) — -, (4.17) 

Vk-i 

k=2 " ^ 

which cancels fl4.15p . The expression (I4.13P is the contribution in the mas- 
ter equation of the particles jumping to the right and leaving the state 
{xn, ■ ■ ■ ,xi) with jump rate a{t)vk, while (I4.14p is the contribution of the 
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particles arriving to the state (xjv, . . . , Xi). Therefore, the jumps to the right 
satisfy the exclusion constraint. 

Secondly, consider all the terms in (14. 6 p which are proportional to b{t). 
They are 



N , N 



^-G(...;t) + ^-G(...,a;, + l,...;t). (4.18) 



Let us denote by m{k) the index of the last particle to the right of particle 
k such that particle m{k) belongs to the same block of particles as particle 
k (we say that two particles are in the same block if between them all sites 
are occupied). Then, (I4.18P takes the form 



N , N 



glHD = - V —G{. . . ; t)+ V —G{. . . , . . . , Xk+k-m{k), . . . ; t). 

•i— ^ 1)7. ^ 1)7 

(4.19) 



tit 



Using (14. lip we get 

—G{. . . , Xfc + 1, Xfc + 1, . . . , Xfe + /c - m{k), ...]t) 

Vk 

= ]_J!l^G{...,Xk + l,Xk + 2,...,Xk + k-m{k),...;t) (4.20) 
= G{. . .,Xk + l,Xk-i + l,. . . ,Xk + k - m{k),. . .]t). (4.21) 

Vk-l 

By iterations we finally obtain 

^1 ^1 
glHD = - V —G{. . . ; t)+ V G'(. . . , Xk+l, Xk-i+l, Xm(k)+l, ...■,t). 

(4.22) 

The first term in (14.220 is the contribution of particles pushing to the left 
and leaving the state {xn, ■ ■ ■ while the second term is the contribution 
of particles arriving at the state {xn, . . . ,Xi) because they were pushed, and 
the particle number k pushes to the left with rate h{t)/vk- □ 

We would like to obtain the joint distribution of particle Nk at time tk for 
A^i > > • • • > A/^m > 1 and < ti < t2 < • • • < tm- By Proposition O, 
this can be written as an appropriate marginal of a product of m determinants 
(by summing over all variables except the Xi''{tk), k = l,...,Nm under 
consideration). 

Notational remark: Below there is an abuse of notation. For example, 
x"(ti) and a;"(ti+i) are considered different variables even if ti = tj+i. One 
could call them simply x"(z) and x"(i + l), but then one loses the connection 
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with the times tj's. In this sense, tj is considered as a symbol, not as a 
number. 

Theorem 4.1. Let us set to = 0, a(to) = b{to) = 0, and Nm+i = 0. The joint 
distribution of PushASEP particles is a marginal of a (generally speaking, 
signed) measure, obtained by summation of the variables in the set 



D = {Xfc(ti), 1 < A; < /, 1 < / < iVi,0 < 2 < m} \ {x^'{ti),l < i < m} 



the range of summation for any variable in this set in Z. Precisely, 
P{xnM = (t^), 1 < i < m\xkiO) =yk,l<k< Ni) 



(4.23) 



const xj^det [^^l_,{x^^ (to))] 



l<k,l<Ni 



D 



X 



n 

i=l 



det[7;^,i,_,(xf'(t,),x^^^(t,_i))]i<fc,K;v, 



Nil 



where 



l<k,l<n 

n=Ni+i+l 



— J) dzz''"^"^e(''(*j)"''(*'^)/'^e(*(*^')"*(*'^)'' 
27ri 

(f) dzz^-y^-\l - vi+^z) ■ ■ • (1 - v^.z), 

2711 



(4.24) 



vl '^l[y>x] and 



n{xl \y) = vl 



(4.25) 

(4.26) 
(4.27) 



The normalizing constant in ^4.24 ) is chosen so that the sum over all vari- 
ables {x^kiti), l<k<l,l<l< Ni, < i < m} equals 1. 

Remark: the variables x""^ participating in the last factor of (14.241) are 
fictitious, cf. (14.271) . and are used for convenience of notation only. 
We illustrate the determinantal structure in Figure [3l 



Proof of Theorem Since the evolution is Markovian, we have 

¥{xnM = Xi\ti), l<i< m|xfc(0) =ykA<k< N^) 
= 5^1(a;t(0)=i/fc,l<fc<iVi) (4.28) 

m 

JJ ¥{xk{t,) = x'liti), l<k< Ni\xk{t^-l) = x5^(t,_i), l<k<Ni) 



X 



i=l 
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t2 



Figure 3: A graphical representation of variables entering in the determi- 
nantal structure, illustrated for m = 2. The wavy lines represents the time 
evolution between and ti and from ti to ^2- The rest is the interlacing struc- 
ture on the variables induced by the det[0„(- ■ ■ )]. The black dots are the only 
variables which are not in the summation set D = D{0)UD*{ti)U- ■ ■UD*{tfn) 
(see Figure m too) . The variables of the border of the interlacing structures 
are explicitly indicated. 



where the sum is over x^{0), 1 < k < Ni, and x'Kti), 1 < k < — 1, 
i = 1, . . . , m. Note that so far the lower index of all variables xf is identically 
equal to 1. 

The continuation of the proof requires a series of Lemmas collected at the 
end of this section, see Section 14.11 We apply Proposition 12.11 to the m + 1 
factors in fl4.28p (including the indicator function, which corresponds to the 
value t = in Proposition 12.11) . Namely, 

F{xk{ti) = x^,{U), l<k< Ni\xk{ti^i) = x'liU.i), l<k<Ni) (4.29) 
= const X ( n ^^^iu)~^-iu-.)'^^ det +^-'(t,) - +^-'=(t,_0, k)] , 

n=l 

where we introduced the notation Oj := a(ti)— a(tj_i), and hi := h{ti) —h{ti^i) . 
First we collect all the factors coming from the We 
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have the factor 



n 



nn 



n / J. J. J. J. 

n=l fc=l n=l 

Vn 



n 

n=l 



n n 

i=l n=Ni+i+l 



Nm 

n=l 



(4.30) 



Then we apply Lemma 14.41 to all the factors det[Ffc ■ ■ )]. For the initial 
condition we have 



^det [F;v,+i-u(2;f(0)- 2/^,0,0 
5(0) 



l<k,l<Ni 



Hdet [Mxt\0),xm)] 



l<k,l<n 



n=2 



(4.31) 



For the other terms, i = 1, . . . , m, we get 



^ det [F^^+i^i^iix^'iti) - x[{ti^i),ai,bi 



D{U) 



X 



Hdet [M4-\U),xnU))] 



l<k,l<n 



(4.32) 



n=2 



Thus, the probability we want to compute in (14.281) is obtained by a 
marginal of a measure on m + 1 interlacing triangles, when we sum over all 
the variables in D{0), D*(ti), . . . , D*(tm), see Figure H] for the definitions of 
these sets. At this point we apply Lemma H75] as follows. For i = 1, . . . , m — 1 
we do the sum over the variables in D{ti). Notice that the remaining variables 
in fl4.30p do not belong to the D{ti), thus we factorize them out. So, r.h.s. 
of (14.281) is, up to a constant, equal to 

Y^m^ X det [F;v,+i-M(xf (0) -i/,,0,0)]^<,^^<^^ 

m-l Ni 
i=0 n=2 

X det ^Fn,+-,+i-i,i{x^"^' (U+i) - x[{ti), ai+i,bi+i) 
J] det [<^„(a;ri(tJ,xr(U)] 



I l<k,l<n 



l<k,l<Ni+i 



X 



I l<fc,Kn 



(4.33) 



n=2 



with the sum is over the variables described just above. By summing over 
the D{ti), the determinant with FAr.^j+i_/^i becomes a determinant with Fi^i 



25 



and the product of the det[v9„(- ■ ■ )] is restricted to n = A^j+i + 1, . . . , Ni. 
Thus, 

g2HD = const x^gSnDxdet [F;vi+i-/,i(xfHO)-yi, 0,0)] 

m 

X n (det (t,) -xf^(t,_i),a„6,)]^<^^^<^^ 

1=1 

X n det [y.„(xr^(t.),xr(t.))] (4.34) 

n=Afi+i+l 

where we set A^m+i = (the contribution from n = 1 is 1). Finally, by using 
Lemma we can include the terms in f l4.30p into the <^n's by modifying the 
last row, i.e., by setting it equal to v^. Thus, 

= const xdet[F^,+i_,,i(xfnO)-2//,0,0)]^^^^,^^^ 

m 

X H (det [Fi,i(xf (t,) -xf'(t,_i),a„6,)]^<^^,<^^ 

n det[0„(a;r^(t,),a;r(t.))]i<,,K„)- (4-35) 



1=1 

X 

n=Ar,+ i+l 



The identification to the expressions in Theorem 14. II uses the representations 
(lOi) and (1331). □ 

The first line represent the initial condition at to = 0, the term with 
in Theorem 14. 1[ These Ni variables evolves until time ti and this is 
represented by the first line (term Tt^ t^). After that, there is a reduction of 
the number of variables from A^^i to N2 by the interlacing structure, which is 
followed by the time evolution from ti to t2. This is repeated m — 1 times. 
Finally it ends with an interlacing structure. If A^i = N2, then the first 
interlacing structure is trivial (not present), while if for example t2 = ti, 
then the time evolution is just the identity. 

In what follows, the picture to keep in mind consists of reading Figure O 
from bottom to top, i.e., in the reversed order with respect to the original 
decomposition. Then nj increases and ti decreases. This corresponds to 
having a sort of vicious walkers with increasing number of walkers when the 
transition is made by the 0's, and with constant number of walkers if the 
transition is the temporal one made by T. 

The measure in fl4.24p is written with the outer product over time mo- 
ments but it can be rewritten by taking the outer product over the index n in 
the variables x'^'s. Let us introduce the following notations. For any level n 
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there is a number c(n) G {0, . . . , m + 1} of products of terms T which are the 
time evolution of n particles between consecutive times in the set {ti, . . . , t^} 
(in other words c{n) is #{i|A^j = n}). Let us denote them by tg < . . . < ^c(„)- 
Notice that = ^(++1), t^' = to, tf^ = ti, and = = t„. Then, the 
measure in (14.241) takes the form 



det[Uxl-\tr'),xnt:^^)))]i<k,i<n (4.36) 



const X Y\ 

n=l 
c(n) 

xndet[7;,,._^(x^(t^),xr(Ci))]i<,,K„ 

o=l 



det[<i_,(xfnO)]i<.,/<^.. 



In Theorem 14.21 we show that a measure on the x^(t^) of the form f l4.36p 
is determinantal and we give the expression for the kernel. Then we partic- 
ularize it in case of the PushASEP with particle dependent jump rates. For 
this purpose, we introduce a couple of notations. For any two time moments 
; ' define the convolution over all the transitions between them by 
0(*ai,ia|) (backwards in time, since forward in the n's). For {ni,ti) -< (^2,^2) 
(see the definition in (12. 5p ). we set 



where 



= t;.! * * r"^+i * ■ ■ ■ * * (4.37) 



T^ = Tt^ t". (4.38) 



For (ni,ti) 7^ {n^i^i) we set 0(*"i'*"2) = 0. Above we used 

^3,t2*^2A =^3A, (4.39) 

which is an immediate corollary of (14.251) . In a more general case considered 
in Theorem 14.2 1 below . if (14.391) does not holds, then T" is just the convolution 
of the transitions between and by definition. Moreover, define the 
matrix M with entries M^^i, 1 < k,l < Ni, 

Mk,i = (0fc * T'^ * . . . * * T^^ * {xl-') (4.40) 

and the vector 

K'5 = <P''''^''*KU- (4.41) 

We remind that the variables x^"^ in M^^i are fictitious, compare with (14.271) . 
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Theorem 4.2. Assume that the matrix M is invertible. Then the normaliz- 
ing constant in ^J^^W is equal to (det M)~^ , the normalized measur^ of the 



form lli4-36 ) viewed as {Ni + . . . + Nm)-point process is determinantal, and 



the correlation kernel can he computed as follows 

/.Til jJlO \ 

K{tllxvXl^2) = -<p'^'^-''-^^\x^,X2) (4.42) 

A'l n2 

+ E E i^m-'Wi * 0(*'(O A^))(xl-\ a:,). 

k=l 1=1 

In the case when the matrix M is upper triangular, there is a simpler way to 
write the kernel. Set 

n 

C5(a;)=E[^~']M(0^*0^*'''''*"O(^^\^) (4.43) 
1=1 

for all n = 1, . . . , Ni and k = 1, . . . ,n. Then, {^^-k} k-i n unique 
basis of the linear span of 

|(0i * 0(*'(i)'*"))(x°,x), . . . , (0„ * 0(*"(")'*"^)(xr\ x)} (4.44) 

that is biorthogonal to 

^<l>"'*"(x)^;'*"(x) =(5„-, i,j = 0,...,n- 1. (4.45) 

The correlation kernel can then be written as 
K{t:i,x,-tZx2) = -0^*"-*"n^i,^2) +EC!5(^i)C3(^2). (4.46) 

k=l 

Moreover, one has the identity 

0^*-'*-^*<:5 = C? (4-47) 
for Hi > n2 and Oi < 02 for rii = 722- 

Proof of Theorem "JTE. The proof is similar to the one of Lemma 3.4 in [4], 



which is in its turn based on the formahsm of [8]. The only place where 
the argument changes substantially is the definition of the matrix L, see [4], 



^With normalized measure we mean that all weights add up to one. If all weights are 
non-negative, it is a probability measure (this is the case for example for PushASEP with 
step initial conditions). 
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formula (3.32). We need to construct the matrix L in such a way that its 
suitable minors reproduce, up to a common constant, the weights (14.361) of 
the measure. Then our measure turns into a conditional L- ensemble in the 
terminology of [8]. 

The variables of interest live in the space 2) = X^^^ U ■■■ UX(^i), with 



("^U-.-UX^g), where xi^^ 



time t" live. Let us also denote / 



= Z is the space where the n variables at 
{1, . . . , A^^i}. Then, the matrix L written 



with the order given by the entries in the set of all variables X 
becomes 



Eo 





El 





E2 


■■ 


Eni-1 





-Ti 











■■ 











-W^[l,2) 








■■ 














-T2 





■■ 

















-^[2,3) 


■■ 

















■ 


-Ts-- 




















■■ 


■ -W^[7Vi-l 

















■■ 





^W) 














•• 






v 

with the matrix blocks in L have the following entries: 

i = n + 1,?/ e X^^+_^^J), 
«G/\{n + l},yGX 












—Tni 




[En]i,y 



(n+1) 
c(n+l)' 



-(n+l) 
c(n+l)' 



and T„ is the matrix made of blocks 



T 



/ Tn,l 













\ 





where 



/U2) 



(4.48) 

(4.49) 
(4.50) 
(4.51) 



X 



GXi"),yGXa 



(4.52) 



(4.53) 



The rest of the proof is along the same lines as that of Lemma 3.4 in [4]. 

Although the argument gives a proof in the case when all variables a;"(tj^) 
vary over finite sets, a simple limiting argument immediately extends the 
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statement to any discrete sets, provided the series that defines / are ab- 
solutely convergent, which is certainly true in our case. □ 

A special case of Theorem 14 .21 is Proposition [HII] stated in Sectional which 
we prove below. 

Proof of Proposition \3.1[ We first prove the statement in the case the jump 
rates are ordered, vi > V2 > . . . , and then use analytic continuation in v/s. 

For Vi > V2 > ■ ■ ■ , the claim is a specialization of Theorem 14. 2 [ The kernel 
depends only on the actual times and particle numbers, therefore we might 
drop the label of t^N Equivalently, we can use the notation (ri,, ti) instead 
of t^'., to go back to the natural notations of the model. For PushASEP we 
have = Fn^+i^i^i{x - yi, 0, 0) and 

Tt^,u{x,y) = - y,a{t,) - a{ti),b{t,) - b{U)). (4.54) 



First of all, we sum over the {xfi(0),l < A; < A^i} variables, smce we are 
not interested in the initial conditions (being fixed). When applied to the 
Fk^i{x, a(ti), b(ti)), the time evolution %.^ti changes it into F^^i^x, a{tj), b{tj)), 

Y,%,,d^.y)FkAy^<ti)MU)) = Fk,i{x,a{t^)Mh))- (4-55) 

This implies that Theorem 14.21 still holds but with t^^ = ti and 

KUi^) = Fn,+i-i,i{x - yi, a(ti), 6(ti)). (4.56) 
We have, see (HSHD, that 

{4>k * Fi^Ni+i-k){x, a, b) = Fi^Ni+2-k{x, a, b). (4.57) 
Using fl4.55p and (14.571) repeatedly one then gets 

= i^7v,+i-^,^,+i-n(a; - m,aitl),bie,)) (4.58) 

which can be rewritten as (13. 4p . 

Next we show that the matrix M is upper triangular. Once again, (I4.55P 
and (14.571) are applied several times, leading to 

Mk,i = 5Z^^Ii^7v.+i-w+i-fc(l/-l/z,«(tcw),K<.)))- (4.59) 

yew. 

Set Ofc = a(tc(fc)) ^k = b{t\k))- Then, 



Mk, = y^vl^j> d..^-^'-V-/-+^-- 1'^'-^^^ . (4.60) 
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(Note that we need the assumption > max{f/};>fc in order for this sum to 
converge.) We divide the sum over y in two regions, {y > 0} and {y < 0}. 
The sum over y > can be taken into the integral provided that \vkz\ < 1 
and then we use JZyyoi'^kzY = jz^- Similarly, the sum over y < is taken 
into the integrals provided that \vkz\ > 1 and one uses ^^y^oivkz] 



y 



For k > I the new term in the denominator, 1 — VkZ, is canceled so that this 
is not a pole and we can deform the contours to be the same. Thus for k > I 
the net result is zero. This is not the case for k < I, since in that case the 
new pole at 1/vk does not have to vanish. The diagonal term is easy to 
compute, since the pole at l/vk is simple. Computing its residue we obtain 
Mk^k = vl'^^e''^''^+''^/''^ and 

det M = JJ y^'+ie^'='^'=+'"=/'"= ^ 0. (4.61) 

k=l 

Next, we need to determine the space Vn where the orthogonalization has 
to be made. We have 



(0, * 0(*^«'*^))(x^\x) = / d^^' 



-x-1 ^ 



(1 - Vk+iz) ■ ■ ■ (1 - vn^z) 
(4.62) 

Once again we divide the sum over {y > 0}, {y < 0} and then deform the 
paths so that the only remaining contribution is the residue at z = 1/vk, 
which is equal to const ■ v^. Thus, Vn = span(t>f , . . . , v^), n = 1, . . . , Ni. 

Finally, we need an expression for the transition between two times, which 
is given by f l4.37p . Every time that we convolute a (pk with T, we get an extra 
factor 1/(1 — Vkz) in the integral. Therefore, if — ^al ^^'^ ^2 > '^i, then 

^«,^4)(x,y) = — i dzz^-y-'^— ^ (4.63) 

while (j)^ ''i' °-2> [x,y) = otherwise. 

Now let us do the analytic continuation. 

As all the functions see (13.41) . can be estimated as 

< const ■ gl^l, xeZ, (4.64) 

for any g > and fi,f2, . • . varying in a compact set, the weights (I4.36P 
can be maj orated by a convergent series for vi,V2, ■ ■ ■ varying in a compact 
set. Further, the normalizing constant (detM)~^ is analytic as long as Vj's 
are nonzero, see (14.611) . Thus, the correlation functions of our measure are 
analytic in vj^s. 
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Set, for = 0, . . . , n — 1, 



fk{x) 



If dz z 



(4.65) 



2ni 7 (1 - Vn_kZ){l - Vn-k+lZ) ■ ■ ■ (1 - VnZ) 



where the integration contour includes the poles v~\, . . . ,v~^ . Note that 
fk{,x) is a linear combination of v^_^, . . . ,v^. Denote by G = [Gk,i]p. n-i 
the Gram matrix 



Since the matrix M is triangular, G is also triangular. Its diagonal elements 
are easy to compute: Gk,k = e°'^^^'"''~^^^^^^^''vl''~^^ . Hence, (14.671) gives a formula 
for $'s that is analytic in f^-'s as long as they stay away from zero. This 
implies that the corresponding expression for the correlation kernel (13. 3p is 
also analytic in f^'s, and thus both sides of the determinantal formula for 
the correlation functions can be analytically continued. Finally, it is not 
difficult to see that the functions (14.651) span the space Ki given by (13.11) . 
which implies the statement. □ 

4.1 Some lemmas 

In this subsection we state and prove the Lemmas used in the proof of The- 
orem 221 

Lemma 4.3. Let us define the function 




(4.66) 



Then for Vi > V2 > 



we have 




(4.67) 




(4.68) 



Then the following recurrence relations holds 



Fk,i+i{x, a, h) = {ipN+i-i * Fk,i){x, a, h) 



(4.69) 



and 



Ffc_i,i(x, a, h) = {ipN+2-k * Fk^i){.x, a, h). 



(4.70) 



From d^. yO| ) and ipn{x, y) = (pn{0, y — x) = {pn{—y, —x) it follows 



Fk-i,i{-x, a, 6) = ^ Fk,i{-y, a, b)ipN+2-k{y, x). 



(4.71) 
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Figure 4: A graphical representation of the summation domains that occurs 
in the next lemmas and theorem. The bold lines passes through the border 
of the domains. 



Proof of Lemma \4-^ We have 

a. = ^ / d^.-VV/- ';-""^>--';-""^"^> . (4.72) 
2vri /r^ (1 - vnz) ■ • • (1 - vm+2-iz) 

Then applying Y.y>x^N'+i-i^'^ = ^^/(l " vn+i-iz) (for \z\ -C 1), we get that 
in the denominator we have an extra factor, which corresponds to increasing 
I by one. Similarly, applying ipM+2-k, the extra factor in the denominator 
cancels the last one in the numerator, thus this is equivalent to decreasing k 
by one. □ 

We define the following domains, which will occurs several times in the 
following. A graphical representation is in Figure HI Let us denote the set of 
interlacing variables at time tj by 

D{U) = {xliU), l<n<Ni,l<k< n\xl+\ti) < x^U) < xl+l{U)}. 

(4.73) 

Then let 

D{U) = {xliU) G D{U)\k > 2}, D{U) = {xliU) e D{U)\n < A^.+i - 1}, 

(4.74) 

and 

D*{U) = D{U) \ {xf (t.)}, D*{U) = D*{U) \ D{U). (4.75) 
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Lemma 4.4. We have the identity 

= const 5^ (n^et [MXk~\ti),^?iti))]i<k,i< 



where we set (fnix^ ,x) = 1 



D{U) "=2 

l<k,l<Ni 

n-1 



X det [Fjv,+i-u(a:fHt.) -xUt.-i), a, ,<^^ (4.76) 



Proof of Lemma 4^ By changing the indices we get that l.h.s. of fl4.76p is, 
up to a sign, equal to 

det [Fjv,+i_i,,(xf -4(^.-i),«,&)]i<fc,KAr, (4-77) 
Using repeatedly the identity fl4.69l) we have 

Fn,kix, a, b) = {LpN,+2-k *---*<fN,,* a, b). (4.78) 

Therefore, 

gTZD = det *---*^N,* F^,+i_;,i)(xf +^-'=(t,) - x[it,^,),a,b)]^^..^^^ 

(4.79) 

We write explicitly the convolution by introducing explicit summation vari- 
ables as follows 

(V5jv,+2-fc *---*ipN,* FN,^+i^i^i){x^'^^'''{ti) - x{{ti^i), a, b) 

= E ( n ^A^.+l-^+n(x:-^-^"(t.), X^t^-'=+"(t.))) 

1<?1<A:-1 

xFN,+i-i,i{x^'{ti) - x[{ti^i),a,b), (4.80) 

where we used the fact that ipm{x, y) = frn{x + c,y + c) for any c E 1i. By 
multi-linearity of the determinant, we can take the sums and the factors ip's 
out of the determinant with the result 

Ni n-1 



x^iU), n=2fc=l 
2<A:<n<Afi 

X det [F;v.+i-u«'(^«) - ^l(^-i)'«'^)]i<,-Kiv, • (4-81) 
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The product of the ip's is non-zero only if ^(tj) < Applying 
Lemma 3.3 in [4] we can further reduce the summation domain to D{ti) 
without changing the result. 

Finally, the product of the determinants of cp^s in the right-hand side of 
fl4.76p is either 1 or depending on whether the variables interlace (belongs 
to D{ti)) or not. This implies (14.761) . □ 



Lemma 4.5. We have the identity 

5^ (Hdet [^„.{xr\t,),xnm 

D{U) "=2 

X det LF7v,+i+i_/,i(a;f'+'(ti+i) -x{{ti),a,h) 



I l<k,l<n 



J l<k,l<Ni+i 



det 



F,^,{xf^'{t,^{)-x';'^^'{ti),a,h) 



(4.82) 



l<k,l<Ni+i 

Proof of Lemma \4.5\ By an analogue (essentially inverse) procedure as in the 
proof of Lemma 14.41 we first get 

Ni+i n-1 

(M) = E {\{\{^n{xt\U).xl^M) 



2<n<Ar,+i-l, 
l<fc<n 



X det 



n=2 k=l 



FAr,+i+i_z,i(a;f'+'(ti+i) - x\{ti), a, h) 



l<k,l<Ni+i 



. (4.83) 



Now we insert by linearity the factor Y[n=i+i Vnixf ^(^i), a;"_,_i(^i)) to terms 
FNi+i+i-i,iix'^''~^^ (ti+i) — x[{ti) , a, h) as well as the sum over the corresponding 
variables. The sums are carried out by using (14. 71 p . from which we get the 



r.h.s. of (021) . 

Lemma 4.6. Let us define 

(f)nix, y) = ifnix, y), (j)nixl-\y) 

Then 



□ 



(4.84) 



<?det [M4-\xr)]-_ 



det [4>n( 



X 



n-l n 



\l<k,l<n ---^ L^"V-fc ^-^l J}l<k,l<n (4.85) 

Proof of Lemma \4.6[ It is a consequence of the fact that both determinants 
are zero if the variables x^ do not interlace and when they do, the matrices 
are upper-triangular with diagonal equal to zero and with equal entries in 
the first n — 1 rows. The only difference is for the last row, where the matrix 
in l.h.s. of fl4.85p has entries 1 and r.h.s. of (14.851) has entries Vn ■ □ 
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5 Asymptotic analysis 



5.1 Flat initial conditions 

To prove Theorem 12.21 we need the uniform convergence of the kernel in 
bounded sets as well as bounds uniform in T. These results are provided in 
the following Propositions 15.11 15.21 15.31 

Let us define the rescaled and conjugate kernel by 

t2{2L+R/2)2X2 

Knm, s,; U2, S2) = K{{m,h),xr; (n^, h), X2)T'/'^^^^^j;^jy^ (5.1) 
where rii = n{ui), ti = t{ui), and 

Xi = [-2ni + Yti-SiT^^% (5.2) 

Proposition 5.1 (Uniform convergence in a bounded set). Fix ui,U2, then 
for any fixed i > 0, the rescaled kernel K^^^ converges uniformly for 
(si,S2)eH,£]2 as 

lim K^r^'"'{ui,si\U2,S2) = S:;^Kj\^^{S^^ui, S;;^si; S^^U2, S:^^S2), (5.3) 

T — >oo 

with Kj^^ the kernel of the Airyi process, see l{2.25\) . and Sy,Sh are defined 
m [Km) . 

Proof of Proposition \5.1\ First we consider the term coming from the second 
integral in (13.231) . namely 

27ri %^ QRt2Z+Lt2/z (^l _ £^ni+n2+xi+l gti (2L+R/2) ' ^^''^^ 

Define the functions 

H{z) = Rz + L/z- {R/2-2L)\n{z), 

go{z) = {n{9) + 9)H{z), 

giiz,u) = -uin'{9) + l)H{z)+u{l-7i'{9))\n{z{l~z)), 

g2{z,u,s) = u\"{9)[H{z)+\n{z{l~z))] + s\n{z), (5.5) 

from which we then set 

foiz) = go{l - z) - go{z), 

fi{z) = gi{l - z,ui) - gi{z,U2) - giil/2,ui) + gi{l/2,U2), 

f2{z) = 5(2(1 Ml, Si) -5(2(2;, M2,S2) -5'2(l/2, Ml, Si) +5(2(1/2, M2,S2), 

fs{z) = -ln(l-z). (5.6) 
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With these notations we get 

The function /o(-z) has a double critical point aX z = 1/2 and the contribution 
for large T will be dominated by the one close z = 1/2. Thus we need to do 
series expansions around the critical point. Computations leads to 

fo{z) = l^,(z-l/2f + 0{{z-l/2Y), 

h{z) = -{ui-u^)K^{z-l/2f + 0{{z-l/2f), 

h{z) = -2{si + S2){z- 1/2) + 0{{z- 1/2)'), 

Uz) = H2) + 0{{z-l/2)) (5.8) 

with 

Ko = 8{8L + R){n{9)+9), = {R + 4L){n'{9) + 1) + 4{1 - n'{e j). (5.9) 

First we choose Fi to be a steep descent pathi for fo{z). Important 
for the later analysis is that the chosen steep descent path is, close to the 
critical point, the steepest descent one. We consider Fi = 7 V 7c V 7, where 
7 = {1/2 + e~™^^^,0 < ^ < 1/2}, 7 its image with respect to complex 
conjugation, and 7c = {1 — l/2e"'^,7r/3 < < 27r — vr/S}. We also have 
fo{z) = Sn{z)R{7ri9) + 9) + Sl{z)L{'k{9) + 9), with 

Sr{z) = l-2^+^ln(^/(l-^)), Sl{z) = ^-l-21n(^/(l-^)). (5.10) 

On 7, simple computations leads to 



dRe(^fl(z)) ^ ^eii+^e) 

de ((i+e) + 2e)((i-0' + 2e)' 

dRe(^L(^)) ^ 64e^((l + 2a^-12a ..... 

de ((i+a+2a'((i-0'+2a' ^ ■ ' 

which are both strictly less than for ^ G (0, 1/2). Now consider the part of 
7c with G [7r/3,7r]. Then 

dRe(gj;(z)) _ 4sin(0)(l-cos(0)) 

d0 5 — 4cos(0) 

dRe(^L(^)) _ 32sin(0)(l-cos(0))(2-cos(0)) 

d0 (5-4cos(0))2 ^ > 



^For an integral / — J^dze*^'^''\ we say that 7 is a steep descent path if (1) Re(/(z)) 
is maximum at some zq £ 7: Re(/(z)) < Re(/(2;o)) for z e 7 \ {zq}, and (2) Re(/(2;)) is 
monotone along 7 except at its maximum point zq and, if 7 is closed, at a point zi where 
the minimum of Re(/) is reached. 
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which are both strictly less than for (f) G [vr/3, vr). For the piece of 7c with 
(f) G (tt, 27r — 7r/3], (15.121) is strictly positive, which is right since when (p 
increases we go closer to the critical point. Therefore the chosen Fi is a steep 
descent path for f'Q^z). 

Take any 5 > and set T{ = {z eTi\\z - 1/2| < 5}. Then, if in (lOl 
we integrate only along F'^ instead of integrating along Fi, the error made is 
just of order 0{e~'^'^) for some c > (more exactly, c ~ (5^ for S small). Thus 
we now consider the integral on F^ only. There, we can use the above series 
expansions to obtain 



.2ri/3 



^^g|Kor(2-l/2)3+(n2-«i)Kir2/3(2-l/2)2-2(si+S2)(2-l/2) 



27ri 

^gO(T(^-l/2)4,T2/3(^-l/2)«,Ti/3(^_i/2)2,(^-l/2)) _ ^^^^^^ 

The difference between (15.131) and the same integral without the error term 
can be bounded by applying — 1| < |x|el^l to 0{---). Thus, this error 
term can be bounded by 

22-1/3 



dz 



g|coKoT(^-l/2)3 + («2-Ul)ciKiT2/3{2-l/2)2-2c2(si+S2)(2-l/2) 



27r 

xO{T{z - l/2)\T^/\z - l/2)^Tl/3(^ - 1/2)^ {z - 1/2)) (5.14) 

for some Co,Ci,C2 which can be taken as close to 1 as needed by setting 6 
small enough. Then, by the change of variable T^^'^(z — 1/2) = w one gets 
that this error term is of order 0{T~^^^) (what is needed is just cq > 0). 

It remains to consider the leading term, namely (I5.13P without the er- 
ror terms. By extending the integral to infinity by continuing the two small 
straight segments forming Ff , the error made is of order 0{e~'^'^). Thus we ob- 
tained that (15.41) is, up to an error 0{e~'^^, T~^/^) uniform for si, S2 G [—i, i]"^, 
equal to 



.2T1/3 
27ri 



^2g3«oT(2-l/2)3+{n2-«i)Kir2/3(2-l/2)2-2(si+S2)(^-l/2)^ (5.15) 



7oo 



where 7oo is a path going from e"^/^oo to e "^/^oo. By the change of variable 
w = {KoTy/^{z - 1/2), we get 

= — [ d«;^e5"''+("2-"i)"''"i/"o^'"'(^i+^2)'"/"o^' (5.16) 



y^Qls~'^{u2-ui)'-^+s:;;^s^\u2-ui)(si+s2) 



38 



with and Sh defined in f l2.13p . Here we used the Airy function represen- 
tation 

— [ di;e"'/3+""'+^'' = Ai(a2-6)exp(2aV3-a6). (5.17) 

To finish the proof, we need to consider the term coming from the first 
integral in 03.231) . namely 

' ( " V^"",^^-^^/-)^^--^^) ^"^"""^^^^^^ (5 18) 
27ri /r/ Vl - (2^+^/2) 2-'i' 

This can be rewritten as 

(|53Sl) = / ^eT2/3(po(»)-po(l/2))+Tl/3(p,(^)_p,(l/2)) 

27ri w; 

with 

Po(w) = {u2-u,){n'{e) + l)H{w)-{u2-Ui){l-n'{e))\n{w{l-w)), 
Pi{w) = -{ul-ul)^^-^[H{w) + ln{w{l-w))]-{s2-s,)ln{w),{5.20) 



where H{w) is the function defined in (15. 5p . Remark that we need to do the 
analysis only for U2> Ui. The function po has critical point at w = 1/2. The 
series expansions of Pq and pi around w = 1/2 are 

Po{w) = p^{l/2) + K^{u2-u^){w-l/2f + 0{{w-l/2f), 
pi{w) = pi{l/2) + 2{si-S2){w-l/2) + 0{{w-l/2f). (5.21) 

We choose as path Fq = {|e"^, G (— vr, vr]}. This is a steep descent path for 
Pq. In fact, for w G Fq, 

Re{H{w)) = (i?/2 + 2L)cos(0) + (i?/2-2L)ln(2), (5.22) 
Re(-ln(w(l - It;))) = ln(2) - In |1 - ti;| = 2 ln(2) - i ln(5 - 4 cos(0)), 

which are decreasing when cos(0) decreases. Thus, we can integrate only on 
Fq = {w G Fo||w — 1/2| < 6} and, for a small S, the error term is just of 
order 0(e~'^^^'^^) with c > (c ~ 5^ as 5 -C 1). The integral over Fq is then 
given by 



-2T1/3 



dwe 



Ki(u2-ui){w-l/2fT^/^+2{si-S2){w-l/2)T'-/3 



2m 
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As above, we use |e^' — 1| < |x|e'^'', to control the difference between fl5.23p 
and the same expression without the error terms. By taking 6 <^ 1 and 
the change of variable w = {w — 1/2)T^/^, we get that this difference is of 
order 0(T~^^^) uniformly for si,S2 in a bounded set. Once we have taken 
away the error terms in (15.231) . we extend the integral to 1/2 ± ioo. By 
this we make only an error of order 0{e~'^^^^^). The integration path can be 
deformed to 1/2 + iIR without passing through any poles, therefore by setting 
w = 1/2 + iyT-^/^ we get 

^ ^ dye-'^^i^---^)y'+'(^^-s2) ^ ^ exp ^ ~ '''' 



^/7TKi{u2 — Ml) V Ki{u2 — Ui 

exp (- ^'''''^'S, ]. (5.24) 



47r(n2 - ni)5, 1 V 4(^2 - mi)^," 

Since all the error terms in the series expansions are uniform for 
(si,S2) G [—iji]"^, the result of the Proposition is proven. □ 

Proposition 5.2 (Bound for the diffusion term of the kernel). 

For any si, S2 G R and U2 — ui > fixed, the bound 



gt2(2L+R/2)2X2 2^1/3 r i 



n2—ni 



< const e^l'i^'^l (^5 25) 

holds for T large enough and consi independent of T . T/ie const is uniform 
in Si, S2 hut not in U2 — Ui. 

Proof of Proposition \5.2[ From the analysis in Proposition 15.11 we just need 
a bound for \s2 — si\ > i, i > fixed. We start with fl5.19p but to obtain a 
decaying bound for large |s2 — si| we consider another path Fq. 

Consider an e with < e <^ 1 and set Tq = {w = pe"^, (f) G [— tt, tt)}, with 

1 + i^pil^, if |s2-Si|<eTV3, 

2 (U2—Ul)Kl ' I ^ -^1 — ' 

-2 + 7=u7, if ^2-^1 > ^rV3, (5.26) 



2 (u2-ui)ki ' 
1 e 



2 («2— 



if S2-S1 < 



We have ^Re{w - iln(u;)) = -psin((/)), ^Re(l/u7 + 21n(u;)) = -^sin(0), 

and ^Re(— \n{w{l — w))) = —jz^^^^y^- Thus Fq is a steep descent path 
for Po{w). Moreover, since on Fq we have Re(ln(w)) = ln(p) is a constant, Fq 
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is also a steep descent path for po{w) plus the term of pi{w) proportional to 
S2 — Si. Let, for a small 6 > fixed, Tq = {w = pe^*^, G (—5, 5)}. Then 

( 157X91) = gr2/3{po(p)-po(i/2))+Ti/3(p,(p)_p^(i/2)) ^^^27) 

\ 27ri u) 

for some c > (for small 5, c ~ 5^). On Fq the Sj-dependent term in 
Re{pi{w) — pi{p)) is equal to zero and the rest is of order Therefore 
the last integral can be bounded by 

mi t ^ -|r2/3(„2-«i)[(7r'(e) + l)(flp+L/p) + (l-7r'(e))p/(l-p)2]</,2+0(;r2/3^4^;rl/3^2) 



-e 

2vr p 

(5.28) 

For 5 small enough, and T large enough, the terms (9(T^/^0^) and (9(T^/^0^) 
are both controlled by the first term in the exponential. Then, by the change 
of variable T^^^cf) = ijj one sees that r.h.s. of (15.281) is bounded by a constant, 
uniformly in T. 

What remains is therefore to bound the first term in the r.h.s. of ( 15.270 . 
By the choice in ( I5.26P of p, \p — 1/2| < e/{{u2 — <^ 1 for e small 

enough which can be still chosen. Series expansion for p close to 1/2 leads to 

Po(p)-Po(l/2) = -2(s2-Si)(p-l/2)Ti/3(l + 0(p-l/2)) 

+ f^i{u2 - ui){p - 1/2)2t2/3(1 + C(p - 1/2)). (5.29) 

By ( 15.261) we obtain the bounds 

Po(p)-Po(l/2) = -l^^^^^{l + 0{e)), if\s2-si\<eT'/M5.30) 
Po(p)-Po(l/2) = - ~ '^'^f '^V + 0(e)), if 1^2 - ^il > eT'/\ 

{U2 - Ui)Ki 



Combining the above result we have 



\mB\< 0(e-=^''') + 0(l)^^ 

(5:31) 

Thus by taking an e small enough and then T large enough the bound (I5.3ip 
implies the statement to be proven, since for any a > 0, there exists a 
Ca< 00 such that e-^^^^-^^)" < C^e-l^^-sii^ □ 



e (i12-"i)ki^ ^ ^' _|_ g (u2-"l)'«l ^ ^ " 
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Proposition 5.3 (Bound on the main term of the kernel). 

Let Ml, U2 be fixed. Then, for any (si, S2) G [—i, 00)'^, the bound 

_2^1/3 r ^Rti(l-z)+Lti/{l-z) ^ni+n2+X2 gf2(2L+_R/2) 2^ 

dz- 



< const e-("i+"2) (5.32) 

holds for T large enough, where const is a constant independent ofT. 

Proof of Proposition \5.3[ Let £ be a constant independent from T, which 
can still be chosen large if needed. For (si,S2) G [— £, the result is a 
consequence of the estimates in the proof of Proposition 15.11 Therefore 
we can consider just (si,S2) G [—£,00)^ \ [— Introduce the notation 
§. = (^Si + i + £)T-2/3, which then belongs to [iP-^/^, 00). 
The integral to be bounded is 



27ri 



£ dze^^°(^)+^'''^i(^)+^'''^^(^)+^=^(^) (5.33) 



where fi{z) and fsi^z) are given in fl5.6p . and fo{z) and f2{z) are just slight 
modifications of the functions in (15. 6p . namely 

fo{z) = {n{e) + e){H{i-z)-H{z)) + hH'^{i-z))-s2H^z), 

f2{z) = g2il-z,uu-i-i)-g2iz,U2,-i-i) (5.34) 
-(72(1/2, m, -£ - £) + (72(1/2, ^2, - i). 

We put si and §2 in fo{z), because they are not restricted to be of order 
T~^/^ (as it was the case in Proposition 15.11) . 

First we need to find a steep descent path for fo{z). We choose it as 
Fi = {1 — pe^"^, G [— TT, tt)} with < p < 1/2, chosen as follows, 

f i-((Si + S2)/fi:o)^/2^ Si + 52<£, /gggx 
^ \ i - (£//^o)^/^ 5i + 52 > e, ^ ' ^ 

for some small £ > to be fixed later. Recall that the Si > 0. 

To see that Fi is a steep descent path, we consider fo{z) term by term. 
First consider G [0, vr], the case (p G [— tt, 0] is obtained by symmetry. The 
term proportional to R{n{9) + 9) satisfies 

-iRed - 2. + 1 !„(./(! = _ P(3- 8pcosW + V) sm(^) ^ „ ,5,3,, 

d0 1 — 2pcos((/)j + p"' 
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for all < p < 1/2, with equality only at = 0, vr. The term proportional 
to L{n{9) + 0) satisfies 

ARe(l/(l-.)-l/.-21u(^/(l-.))) = J{^-M) + 2p^f-£)M^) , 

d(j) ^ ' ' ^ (1 -2pcos(0) +p2)2p 

(5.37) 

for all < p < 1/2, with equality only at = 0,7r. Finally, Re(ln(l — z)) 
is constant on Fi and — Re(ln(2z)) = — ln(2|2;|) is strictly decreasing while 
moving on Fi with |0| increasing. 

For a small 5 > 0, Ff = {1 - pe'"^, G (-5, 5)}. We also define 

g(p) = exp (Re(T/o(l - p) + T^/^ill - p) + T^'^ - P))) ■ (5.38) 
Since Fi is a steep descent path of /o(-2), the integral over Fi \ Ff is bounded 

by 

g(p)0(e-'=^) (5.39) 

for some c > independent of T. The contribution of the integral over V\ is 
bounded by 



Q{p) 



27ri 



d^e2^{/0(^)-/o(l-p))+T2/3(/i(^)-/l(l™p))+Tl/3(/2(^)_/2(l™p))+/3(^) 

(5.40) 



The series expansion around <^ = is 

Re(/o(l - pe^^) - /o(l - p)) = -7i0'(l + 0(0)) (5.41) 

with 



g^p (7r(e) 
7i = izr. ^ ^ 



+ ^)(l-2p) / i?p(3-2p) L(l-p + 2p^) \ 

;i-p)2 V 4 ^ 3p(i-p) ; 



2(1 -p)2 

(5.42) 

and 

Re(/i(l - pe'<^) - /i(l - p)) = 720'(1 + O(0)), (5.43) 

with 

72 = (n2-ni)/«i + 0(p- 1/2). (5.44) 

Finally, Re(/2(1 — pe"^) — /2(1 — p)) = C(0^). Thus, by the change of variable 
z = \ — pe^"^, the above estimates, and by setting 7 = 71 + 72T~^/^, we get 

(I530D = g(p)^^ d0e-^'^'^(^+'^(^))(^+^(^'^'')). (5.45) 
2vr(l - p) 



43 



By choosing 6 small enough (independent of T) and then T large enough, 
the factors with the error terms can be replaced by 1/2, thus 

Remark that, the worse case is when 7 becomes small, and this happens 
when p — > 1/2, i.e., it is the case of small values of Si + 52- But even in 
this case, jiT^/^ ~ (si + S2 + 2i + 2ifl'^ > {Ufl'^. Since 72 is of order one, 
^2^1/3 _ 2^1/3 _|_ > for £ large enough. So, for by setting i large enough, 
f l5.40p < const (5(p). This estimate, combined with (15.391) . implies that the 
Proposition will be proven by showing that (5(p) < const e~^*^"*"'^^\ Since 
1 — p is close to 1/2, we can apply the series expansion of /j around z = 1/2. 
The expansion of /i is in (15. 8p . while the one of /2 is the same as in (15. 8p 
with si + S2 = -2i - 2i. Finally, 

/o(z) = Iko{z - 1/2)3(1 + 0{z - 1/2)2) _ ^ _ ^/2)(i + 0{z - 1/2)). 

(5.47) 

First consider Si + S2 < s. Then, with p chosen as in (I5.35p . we get 

^g-2{i;+£){Si+S2)/.-i/W3(i+o(^)) 

= e~i(*i+^2+2£+2l)3/2K-i/2(i+o(^))^(u2_„j)Ki(si+s2+2£+2l)K-i{l+0(v^)) 

^g-2(^+Z)(si+.2+2^+2£>-^''"T-V3(i4.ci(Vi))_ ^5_4g^) 

Recall that si + S2 + 2i+2£ > 2£ > 1 for £ > 1. Therefore by choosing i large 
enough (depending only on the coefficients kq, ki,ui,U2 which are however 
fixed), all the terms are controlled by the first one, i.e., 

Qip) < e-h'i+s2+2i+2ef^^.-'/' < ^_i(,,+,,)3/2,-i/2^ ^^^^^^ 

Since this decays more rapidly that exp(— (si + S2)), the Proposition holds 
for Si + S2 < e. 

The last case is Si + §2 > e. In this case, with p chosen as in (15.350 . we 
obtain 

^g--4{f+f)«,7^/2eTi/3(i+ci(^))^ ^^^^^^ 

But now, e/3 — {§1 + §2) < ~f (^1 + ^2), thus the first term in the exponential 
is, up to a positive constant,— A/eT-^/3(si + S2 + 2£ + 2i), which dominates 
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the second term ~ eT"^^'^ < Si + S2 + 2£ + 2£, and it also dominates the third 
term. Therefore, for any choice of e and £ made before, we can take T large 
enough such that 

Q(p) <e-^^^''"(^^+^2), (5.51) 
which ends the proof of the Proposition. □ 

Proof of Theorem \2.S[ The proof of Theorem 12.21 is the complete analogue 
of Theorem 2.5 in [3]. The results in Propositions 5.1,5.3,5.4, and 5.5 in [3] 
are replaced by the ones in Proposition 15.11 15.21 15.31 The strategy is to write 
the Fredholm series of the expression for finite T and, by using the bounds 
in Propositions 15.21 and 15. 3[ see that it is bounded by a T-independent and 
integrable function. Once this is proven, one can exchange the sums/integrals 
and the T —>■ 00 limit by the theorem of dominated convergence. For details, 
see Theorem 2.5 in [3]. □ 



5.2 Sketch of the result ([233]) 



With the rescaling (12.91) and (I2.19p . the rescaled kernel writes 

K'^'^im, su U2, S2) = K{{mM), xi; (ris, t2),X2)T'/\ (5.52) 
The main part of the kernel (the second term in (13.171) ) writes 

7-1/3 f r gr/o(t«)+T2/3/i(M,;«i)+Tl/3/2(^„;„j,Si) ^ 



W Ui 



(27ri)2 Jr^ Jp^ e^^o(^)+'^''''^i(^'«2)+r''''/2(^;«2,s2) w{w - z) 

with 

fo{w) = [nid) + 6) [Rw + ^) + {n{d) - 6) In (1^) - a, In(^), 
f^{w- u,) = - [(7r'(^^) + 1) [Rw + ^) + (7r'(^) - 1) In (i^) - a, ln(. 
f2{w- M„ Si) = [\iT"{e) + ^ + In (i^)) - aa] uj + ln(w;). (5.54) 

The parameter /i is actually the position of the double critical point of 
fo{w). Series expansions gives 

/o(^) = fM-^{w-^^f + o{{w-^^Y), 

fi{w]Ui) = fi{n;ui) - uiKi{w - fif + 0{{w - fif), (5.55) 
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The terms /i(yu; Ui) and /2(/x; Ui, Si) cancel out by an appropriate conjugation 
of the kernel (I5.53I) . We denote by ~ an equality up to conjugation. Thus, 
asymptotically, (15.531) goes to 

rpl/3 I" dz e^i^o{w-lJ.)^T/3-uiKi{w-ij.y^T^/^+T^/^{w-ij.){si/ti-Klul/Ko) 

dw 



/i(27ri)2 J^^w- Z e-«^o(2-/^)»T/3-«2«;i(^-M)2T2/3+Ti/3(2_^)(s2//,_«:2„2/^^) 

(5.56) 

With the change of variable {w — fi)(KQTY^^ = W, (z — jj){KQTY/'^ = Z, we 
then obtain 



W - Z ^^W3+UlW^Kl/Kl^^-W{si/tl-Klul/K0)/Kl/^ ' 

(5.57) 

Let us denote by 5*1, = fiK^^^ and Sh = Hi^n^^^ the vertical and horizontal 
scaling. Then 

(ESTD = S;'KA,{Sh'ui, S;hi; S^'u2, S;\s2) (5.58) 

where Kj^^ is the extended Airy kernel associated to the Airy2 process. An 
asymptotic analysis of large deviations similar to Propositions 15.21 and 15.31 
above would then lead to the result of (12.231) . 
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